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ABSTRACT 


The effects of radiation on the performance of modern rocket pro- 
pulsion systems operating at high pressure and temperature have been 
recognized as a key issue in the design and operation of various liquid 
rocket engines of the current and future generations. 

Critical problem areas of radiation coupled with combustion of 
bipropellants are assessed and accounted for in the formulation of a 
universal scaling law incorporated with a radiation-enhanced vaporiza- 
tion combustion model. Numerical algorithms are developed and the per- 
taining data of the Variable Thrust Engine (VTE) and Space Shuttle Main 
Engine (SSME) are used to conduct parametric sensitivity studies to 
predict the principal intercoupling effects of radiation. The analysis 
reveals that low enthalpy engines, such as the VTE, are vulnerable to a 
substantial performance set back by the radiative loss, whereas the 
performance of high enthalpy engines such as the SSME, are hardly 
affected over a broad range of engine operation. Additionally, combus- 
tion enhancement by the radiative heating of the propellant has a 
significant impact in those propellants with high absorptivity. Finally, 
the areas of research related with radiation phenomena in bipropellant 
engines are identified. 
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1 . INTRODUCTION 


The performance of modern rocket engines (Sutton and Ross, 1975) 
is primarily affected by the competitive rate processes involving multi- 
phase turbulent reacting flow with phase change and heat transfer at an 
elevated chamber pressure, typically at 10 3 psia and chamber tempera- 
tures of 5500 to 6500°R at which levels the flame emission (Afgan et al. 
1974) due to the molecular gas radiation plays two significant roles. 
Firstly, the radiation emission leaving the system constitutes a thermal 
energy loss and, thus, directly and indirectly affects the major per- 
formance parameters through the chamber temperature. These performance 
parameters include critical mass flux, specific impulse C*, and total 
thrust. Secondly, the radiative heating of propellants enhances the 
vaporization and combustion of droplets and, thus, the process serves 
to improve the overall combustion efficiency of the motors. 

Because of these two diametrically opposite roles played by radia- 
tion, the overall impact of flame emission may result in performance 
loss or gain depending on the relative preponderance of heat loss rela- 
tive to propellant gasification rate enhancement. This constitutes the 
focal point of quantitative scrutinization. 

The principal factors affecting the nature and extent of the 
impacts of radiation are traced down to the fundamental physics of 
emission, absorption, and scattering of gas molecules, droplets, soot 
particles in hydrocarbon propellants, and the combined radiation and 
conduction-convection heat transfer in the chamber (Summarf ield, 1960). 

The practical significance of the above-mentioned factors depends on 
the basic design and operating parameters including propellant properties, 
mass flow rates, mixture ratio, chamber pressure, spray characteristics, 
and the combustor size. 

One of the critical design issues of the current and future rocket 
engines is the lack of a comprehensive criterion that predicts and 
classifies the motor performance characteristics, in the presence of 
the radiation processes, in terms of the design and operating parameters 
described above. 

Review of the literature accumulated over the past several decades 
reveals that there has been little or no in-depth research which pertains 
to the fundamental issues of radiation in liquid rocket combustion 
chambers. Furthermore, there has been a genuine lack of vital engineer- 
ing data, modern analytical tools, and computational fluid dynamics 
(CFD) required for the analysis of the bipropellant and monopropellant 
combustion processes. In fact, many of the existing data and analytical 
tools developed for the radiative treatment of large scale industrial 
combustion appliances (Afgan et al., 1974) and air breathing engines 
are incomplete and inadequate for application in rocket motors for 
performance and design optimization. To this end, a comprehensive 
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motor performance scaling law is developed, on the basis of rudimental 
thermodynamic and combustion principles, to aid in assessment and iden- 
tification of the basic issues and qualitative prediction of the motor 
performance characteristics of modern liquid rocket engines. 

The objectives of this study are (1) to develop an algebraic 
scaling law of motor performance and propellant droplet gasification 
laws, (2) to assess and compare the performance of the VTE and SSME with 
respect to regenerative heat recovery systems in each type of engine, 
and (3) to identify specific technical problem areas where radiation 
modeling and numerical simulation of radiation-coupled combustion 
processes are required. 
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2. SCALING LAWS OF ROCKET MOTORS 


A performance scaling law is developed with (1) thermodynamic 
relationships that determine the chamber temperature as a function of 
propellant properties, chamber inlet conditions, wall heat transfer by 
conduction, convection, and radiation, combustion and vaporization 
efficiencies, and a regenerative efficiency; and (2) propellant droplet 
laws that predict the rate of vaporization and combustion under the 
effects of combined radiation and conductive heat transfer. 

Such a law must preserve the universality for any type of propellants 
and engine configuration, regenerative or non-regenerative, and all 
engine power levels. 


2.1 Thermodynamic Performance Model 


A schematic. Figure 1 (a and b) , shows a regenerative rocket engine 

that admits fuel and oxidizer at flow rate nw, and m n , respectively, at 

r U 

an effective inlet temperature T^. The reheater receives heat at the 

rate of q_, radiative heat, and q r , conductive heat, from the bulk of 

the gas in the combustion chamber. The fraction of heat received, n n 

(q_ + q ) , is used in propellant preheating and the remaining part is 
K L 

lost to the environment. The bipropellants enter the motor at a tempera- 
ture T^ 1 and gasify at the efficiency of n RVR for fuel, and n 0UR for 

oxidizer, respectively. The overall combustion efficiency is The 


subscript R reflects the inclusion of the radiation effects during the 
phase change and propellant combustion. The heat released by combustion 
raises the temperature of the gas product, T^, and partly transfers heat 

to the wall by conduction q c and radiation q R . Thus, T^ is, in general, 

lower than the adiabatic flame temperature. The flame emission from the 
hot combustion zone is assumed to be uniformly distributed in the 
chamber and is transported to the wall, held at the temperature T^. 


2.1.1 Scaling Law of Rocket Motor Performance 

The overall energy balance of a rocket motor is given by 


A 


f (l+<t>) C p (T c -T.') « -q c - q R + n c A f (!+<)>) Q c 


m f ^VFR L F + ^ n VOR L 0- ) 


( 1 ) 
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Fig. 1(a). Schematic of Liquid Propulsion System 
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Fig. 1(b). Schematic of Thrust Chamber 
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and for the reheater 

\ (q C + q R } = A f (1+<f)) C P (T i' " V (2) 

where the wall heat transfer rates q and q are 

C R 


71 N k L 

q c V— < T C - V < 3 > 

1r ' 0 ’ DL (e g I c 4 - % T w 4 ) (4) 


(J) is the mixture ratio, N u is the Nusselt number, k is the thermal con- 
ductivity of the gas, e is the emissivity, a is the Stefan-Boltzmann 
constant, a is the surface absorptivity, D is the diameter, and L is 
the length of the chamber, respectively. 

By eliminating ' between equations (1) and (2) and by adopting 

a non-dimensionalization of the resulting equation with appropriate 
reference properties, one obtains the following algebraic scaling law: 


* e c 


4 


+ 9, 


= K 


where 

e (l - n h ) 

Y = 1 + a (1 - q h ) 

K - 1 + 9 i + 01 (1 ~ 9 w + $ W (1 - V 9 w 4 

1 + a (1 - n h ) 

R k it L N u k ir L C* 

2 C p m f (1+*) 2 C P P C A t 


(5) 


( 6 ) 

(7) 

( 8 ) 


V 1 1 1 - 5 


(9a) 


3 3 

aeirDLQ., a e it D L Q„ C* 

£j _ JCj 

C p 4 m f (1-N.) 4 Cp 4 P c A t (1-H>) 3 




a ay ir D L Q e 3 c ^ i D L Q £ 3 C* 

Cp 4 m f (l+4>) 4 Cp 4 P c A fc ( l+4>) 3 


(9b) 


e 


(1-H>) C p T 
% 


( 10 ) 


Qg (1-N>) n c Q c ^VFR L F + ^ n V0R L (P 

in which Cp is the specific heat of the combustion product, C* is the 

characteristic exhaust velocity, A is the throat area, L_ and L-. are the 

t r u 

latent heat of vaporization of the fuel and oxidizer, respectively, (j> is 
the mixture ratio, P^, is the chamber pressure, and k is the thermal con- 
ductivity of the gas. 


2.1.2 Conductive and Convective Heat Transfer 


The conductive and convective heat transfer in rocket motors has 
been correlat ed by various semi-empirical expressions involving the 
Nusselt number. For example, the Colburn equation gives an expression 
of the Nussej t number in terms of the Reynolds number and Prandtl 
number. In some sophisticated empirical laws, the effects of the 
boundary layer temperature gradient on the various gas properties near 
the chamber wall are accounted for. For example, Bartz (1958) gives the 
following semi-empirical formula: 


Nu 




( 12 ) 


The subscript m refers to properties evaluated at the arithmetic mean 
temperature of the local free-stream static temperature and wall tem- 
perature, pg is the local gas density of the free stream, v is the local 
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gas velocity, subscript c refers to mean chamber temperature, and k is 
the thermal conductivity of the gas. 

In the present analysis, all the local values will be approximated 
by appropriate mean values. 


2.1.3 Radiative Heat Transfer 

The radiative emission in a combustion chamber is largely due to 
molecular rotational and interatomic vibrations of heteropolar gaseous 
molecules such as water vapor, carbon monoxide, carbon dioxide, hydro- 
carbons, ammonia, and nitrogen oxides. Symmetrical molecules such as 
hydrogen, oxygen, and nitrogen show no appreciable emission bands in 
those wave regions important for chamber radiant heat transfer. In a 
heterogeneous mixture, such as bipropellant sprays, clouds of finely- 
spread disperse phases, such as soot particles of hydrocarbon-based 
propellant combustion, the radiation power could increase by a large 
factor of 2 to 10 times the equivalent homogeneous gas mixture. 

The radiant heat transfer to the wall by the hot gas emission is 
given by 


Or - 0 4 (E g t4 - “w t w 4) (13) 

where e is Che overall gas emissivity, and a w is the effective surface 
absorptivity of the wall. If the wall surface temperature is close 
to the gas temperature T, then a w approaches Eg. When and T are 
not too far apart, the following approximation may be used: 

e g = ay " e gav (4 + a + b - c)/4 , (14) 

where e is the gas emissivity evaluated at the arithmetic mean 

g 3 av 

temperature of T and T^, and the coefficients a, b, and c are given by 


3 (In e g ) 

p g V 


(15) 


3 (In Eg) 
8 ( Ln T) 


(16) 
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c 


(17) 


9 (In 3 ^) 3 (In e g ) 

= 9 (In T) + 9 (In T) 


where L R is the radiation path length of the flame and P is the gas 

1'" O 

pressure. Numerical data of the coefficients a, b, and c for CC^ and 
^0 have been published by Hottel (1959). 


In the present analysis, the numerical value of absorptivity, a 


w 

at the wall that appears in equation (4) can be taken to be different 

from the gas emissivity value e q . This difference is incorporated in 

the definitions of the expression of g which contains e and of g which 

w 

contains a . Further assumptions concerning the numerical values of e 
w g 

and a are discussed in sections 3 and 4. 
w 


2.1.4 Performance Profile of the Combustion Process 

The method of application of the scaling law, equation (5), for the 
evaluation of combustor performance is described below. 

2. 1.4.1 Non-Adiabatic and Adiabatic Chamber Temperature 

The non-dimensional chamber temperature 0^ is calculated from the 

scaling law for the prescribed engine design and operating parameters. 

By using equations (8) to (18), one can, in general, determine the 

values of a, g, g , and Q to estimate the numerical values of K and y 
w h 

from equations (6) and (7). The dimensionless temperature 0 is then 

calculated from equation (7) and plotted as a function of y and K (Fig. 
2). The actual temperature T„ is finally determined from equations 
(10) and (11). 

The non-dimensional isothermal lines, i.e. constant value of 0^,, 

are a family of straight lines that intersect the K-coordinate at the 
chamber temperature 8^, (Fig. 2). Note that when = 1.0, the non- 

dimensional temperature 0^ assumes the value of 1 + 0^ for any value of 

a and g. This value is equal to the "adiabatic flame temperature." 

An examination of the asymptotic behavior of the scaling law at 
larger and smaller values of y reveals a rather interesting 
feature of the non-adiabatic flame temperature profiles. When y is 
much smaller than unity, and K is close to a constant value, i.e., a, 
g, and 3 W << 1 5 the variation in the chamber temperature with respect 


VIII -8 



ORIGINAL PAGE IS 
OF POOR QUALITY 



VIII -9 




ORIGINAL PAGE IS 
OF POOR QUALITY 


to the variations in the numerical values of y* <*, and 3 W is negli- 
gibly small. This behavior is similar to what may be termed an Asymp- 
totically adiabatic" pattern wherein the chamber temperature is close 
to the adiabatic flame temperature in the broad range of the values of 
a, 3, and 3 • Physically, the case corresponds to the rocket engines 

with a "high enthalpy throughput" in comparison to the combined radiative, 
conductive, and convective loss, so that the latter thermal energy loss 
hardly affects the chamber temperature. On the contrary, when the values 
of y are not small, the non-dimensional chamber temperature variation, 
associated with the variations in the values of a, 3, and 3 W > is not 

negligibly small, i.e., the chamber temperature depends strongly on the 
loss mechanisms. This corresponds to those "low enthalpy throughput" 
engines in which the combined heat loss by radiation, conduction, and 
convection constitutes a significant fraction of the enthalpy flow rate 
through the engines. 


2.1.5 Engine Classification, High and Low Enthalpy Engines 

Two remarkably different chamber temperature characteristics at low 
Y values (y « 1) corresponding to high enthalpy flow and high y (y -1 1) 
for low enthalpy flow, suggest unique classification of all the liquid 
rockets into high and low enthalpy engines. In order to illustrate 
this feature. Table 1 gives a set of data of the SSME and VTE used to 
estimate the values of y of each engine. 

TABLE 1. VTE and SSME Data 




VTE 

SSME 

Combustor 

Diameter (ft) 

0. 

204 

1 

.40 


Length (ft) 

0.333 

0 

.67 

Propellants 

Propellant 

Monomethyl- 

Hydrazine 

Nitrogen 

Tetroxide 

Hydrogen 

Oxygen 


State at inlet 

Liquid 

Liquid 

Gas 

Liquid 


Boiling Point (°R) 

649 

530 

36.9 

163 


Heat of vapori nation 
(Btu/lb) 

377 

178 

190 

85 


Specific heat 
(Btu/lb °R) 

0.7 

0.378 

3.75 

0.5 


Mass flow rate 
(lb/sec) 

0.177 

0.2932 

150 

900 


Inlet temperature 
( e R) 

460 

460 

500 

163 


Heat of combustion 
(Btu/lb) 

8910 

112 

!50 

Scaling 

Parameter 

Y 

-3 

10 J 'L 

1* 

10' 5 -V 

10- 3 * 


K 

1.02 'v 

1.07 

1.03 a. 

1.04 


* The upper value is estimated at N - 300, and the lower value is estimated at N - 3 

u 
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In Table 1, the values of y for the SSME are much smaller than those 


-2 -3 

of the VTE by an order of magnitude of 10 to 10 . This is primarily 
caused by the significant difference in the total mass flow rate of the 
SSME, i.e., mp + * o = 1050 lb/sec, and the flow rate for the VTE which 

is 0.470 lb/sec. The latter value amounts to only 0.047 percent of the 
SSME. 


The consequence of the difference in y values is the dramatic 
difference in the chamber temperature. As a rule of thumb, the chamber 

_3 

temperature reduction in high enthalpy engines with y £k 10 , is of the 

order of "a degree," and that in low enthalpy engines with y = 10 ^ is 

2 3 

of the order of 'vlO to 10 degrees. The detailed chamber temperature 
loss in the SSME and the VTE is discussed at length in Section 4. 


VIII-11 



3 . 


PROPELLANT DROPLET VAPORIZATION AND COMBUSTION IN THE RADIATION 
ENVIRONMENT OF A COMBUSTION CHAMBER 


The effects of the radiation enhancement on the propellant vaporiza- 
tion, combustion, overall combustion efficiency, and, subsequently, the 
chamber temperature, are included in the scaling law using the radiation 
modulated efficiencies n^ R » and n Q y R . However, the values of these 


efficiencies can presently only be estimated because of the lack of 
analytical or experimental data and formulas that interrelate the values 
of these efficiencies with a combustor design and operating condition. 


The basic steps required in the determination of these efficiencies 
are (1) the acquisition of the droplet laws that account for the radia- 
tion enhanced vaporization and combustion, and (2) the prediction of the 
detailed spray combustion field by a comprehensive CFD code that incorp- 
orates the radiation modulated droplet laws. In this section, the first 
step extends the classical single droplet theory (Williams, 1985) and the 
second step describes an ad-hoc physical model that interrelates the 
overall combustion and vaporization efficiencies with the single droplet 
laws obtained in step 1 . 

The proposed model of radiation-enhanced droplet vaporization and 
combustion laws is developed on the basis of the classical single droplet 
theory with the following additional features. The radiation environment 
is isotropic and the gas is assumed to be a gray medium. The droplet has 
a constant absorptivity and the radiative transfer is in radiative non- 
equilibrium. The state of non-equilibrium radiative heat transfer will 
therefore require the treatment of the combined radiation-conduction 
convection mechanism via the solution of radiative transfer equation 
coupled with the conservation laws of mass, momentum, and energy of the 
participating gas flow field. 


3.1 Mathematical Formulation 

The conservation laws governing the droplet processes (Fig. 3) are 
given by 

"T IF (r2 V = B k W k = m ’ °» F ’ T (18) 

r 

where 

■ p v j = F Fuel (19) 
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Fig. 3. Single Droplet Combustion in Radiative Environment 


da . 

= K a j " PD 


j = o oxidizer 


( 20 ) 


*T = a T 



+ q R 


( 21 ) 


in which p is the gas density, v is the gas mixture velocity, D and X 
are mass diffusivity and thermal conductivity of the gas mixture, respec- 
tively, Cp is the gas specific heat at constant pressure, and q R is the 

radiative heat flux, a. and a_ are Schvab-Zeldovich variables defined by 

j T 7 



"-v.’) W. 
J J 


a 


T 



qo 


( 22 ) 


(23) 


q R 



(24) 


where y ^ is the mass concentration, W. is the molecular weight of the jth 
gaseous species, and v." and v.' are the stoichiometric coefficients 

J J 

which appear in the one-step reaction process described as follows: 


U5) 


and q° is the heat of combustion 

q° =£h.° W. (V i * - v ± ") (26) 

in which h/ 5 is the standard heat of formation per unit mass for species 
i at reference temperature. Q R and q R are the dimensional and non- 
dimensional heat fluxes, respectively. 
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The radiative heat flux is given in terms of the radiation 
intensity 1^ as follows 


n 1 

J o J -l 


q r = 2ir 


VI I x (t x , y) dy dX 


( 27 ) 


and 1^ is governed by the following radiative transfer equation 


91, (lV) 91, 

y -5 — + ; 5 — + I, = (1-w. ) I, , + 

9x^ xX 9y X X' bX 


\ r 1 

2 l-i 


I A (T A ,y) dfi 


(28) 


Here is the optical thickness, ft is the solid angle extended by a 

surface relative to an observer, X is the wave length, y is the cos 0 
with 0 being equal to the angle extended by the position vector in radial 
direction, and that of the radiation beam. The terms a^, <5^, and are 

absorption, scattering, and extinction coefficients, respectively. 8^ 
and w^ are given by the following expressions 



w 


X 


V s 


A 


The subscript b appearing in the radiation intensity refers to that 
of black body radiation. 


The boundary conditions of the present problem are (1) an impermeable 
droplet surface with respect to oxidizer and combustion product, ( 2) 
the balance of the combined radiation and conductive heat transfer with 
the latent heat of vaporizing species leaving the droplet, (3) inten- 
sities of radiation on the droplet surface and at the gaseous environment 
at r oo, and (4) the temperature and concentration of species at r 
These conditions are summarized in the following: 


r 




VI 1 1 — 15 


4ir 


2 (dT n \ 

r * fa • U. ■ " 


n 2 „ T ' 4 

1 <V> = —5 


r + « : 


T - T » . Y i - Y i» 


I (®,y) = 


2 n T 4 
n a T 


(30) 


(31) 


(32) 


(33) 


3.2 Radiation-Enhanced Vaporization, Burning Laws and Flame Location 

By following the mathematical procedure adopted in the classical 
single droplet theory, one can predict the vaporization or burning rate 
as the characteristic values of the two-point boundary value problems. 

The results of the analysis are listed in the following: 

(1) For a vaporizing droplet, the vaporization rate &y R is calcu- 
lated in terms of the "radiation-enhanced vaporization transfer number," 
By R , as follows 


&VR = 4 ir pD r^ In (1 + B vr ) 

.T 


/, 


C P dT £ 

T t . .0 eH 


b vr l 


m q 


4tt (pD) L J o 


L 


q R (5’K' 2 e 5 ' d5' 


where 


= m f 

•'n 


dr 


o 4ir pD r 


2 


(34) 


(35) 


(36) 


(2) For a combusting droplet, the burning rate and the flame loca- 
tion r^ are predicted to be 
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m CR = 4 tt pD In (1 + B CR ) 


(37) 


and 



In 


In (1 + B cr ) 


(■* 


W F V F 

W v 
o o 



(38) 


where the radiation-enhanced combustion transfer number, 
expressed by 


B CR’ 


is 


CR 



o 

m q 

4it (pD)^ 


/ q R (S') V 2 e ? ' d? 


j L. (39) 


Note that when the temperature of the environment is higher than 
that of the droplet surface, the radiative heat flux points toward the 
droplet, thus, the q D has a negative value. Under such circumstances, 

equations (35) and (39) show that the transfer numbers B VR and B CR are 

larger than the values determined for the cases when the radiation is 
absent. 

2 

A special case in which q R r = constant, equations (34) and (35), 
reduce to the following form: 

*VR or CR " 4* pD r £ In (1 + or CR ) (40 


where 


VR 


1 - 


(Q R */mL) 


and 


"CR 


1 - 


(Q R */mL) 


(41) 


(42) 
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Here Q R * is the rate of the radiative heat transfer to the droplet 
surface, i.e.. 


V ■ 4 * x l q r Cr ' r i> 


(43) 


Thus, the enhancement of the vaporization or combustion of a droplet can 
be expressed by the transfer number factor given by 


_ b vr _ b cr 1 

R B c 1 - (Q R */mL) 


(44) 


where e R is greater than unity. 

An alternative expression that describes the vaporization and com- 
bustion enhancement is to use the radiation factor defined by 



(45) 


The vaporization rate can be expressed in terms of the radiation 
factor n R as 


= 4ir pD r £ In 1 



(46) 


The droplet vaporization and combustion radiation enhancement factors 
Cy R and C^, R can be defined as the ratio of the vaporization and combustion 

rates as follows: 


C 


VR 


*VR 

*V 


(47) 
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(48) 


m. 


CR 


"CR m. 


where lii^ and are the vaporization and combustion rate for the same 

environment conditions at r + «, but in the absence of radiative effects 
on the droplet rate processes. 


3.3 Radiation Flux and Intensity 

The dimensionless radiation heat flux q R appearing in the vaporiza- 
tion and combustion rates, equations (34), (35), (37), and (39), is a 
function of the intensity of radiation 1^, equation (27). Thus, the 

determination of the heat flux requires the knowledge of the solution 
of equation (28). 

According to the assumption of a gray gas and an isotropic scatter- 
ing, the solution of the radiative transfer equation can be expressed in 
the following power series solution involving Legendre functions 


I(t,w) 


T 2m+1 p 

4ir m 

m=0 


(y) (t) 


(49) 


where y = cos 0. By substituting equation (49) into equation (28), 
which is specialized to the gray gas approximation, and by equating the 
coefficient appearing in each order of the Legendre function, one obtains 
the following set of equations: 


m = 0 ; -j^r + 2 — + (1-w) ip Q = 4ir (1-w) I b ~ j w 


(50) 


dip. 


'Po dip. 


m = 1 ; 2 -r^ +6—4- + 3 jfi, =0 

dr t dr r l 


(51) 


m = n ; (n+1) 


[V- 


(n+2) 


*n+l 


+ (2n+l) tj> 


n 


F^n-l , .. *n-ll 

n [-jp- - (n-l) — J - 


(52) 
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The net radiation flux q R * and incident radiation G can be expressed 
in terms of the integral as follows: 



U I(t,p) du s 


I(t,y) du s ^ 


1 


(53) 


(54) 


where the first two terms ij; o and appearing in equation (49) are chosen 
to construct the heat flux q R * and incident radiation G*. 

Upon substituting equations (53) and (54) , one obtains the following 
equations governing q R *: 


d2 V . 2 d q R* 

d, 2 * d ^ 


3 (1-w) 




4ir 


(1-w) 


lib 

dT 


(55) 


and 



1 d G* 
3 dt 


(56) 


In general, equation (55) must be solved in conjunction with the 
conservation equations (18) together with the boundary conditions, equa- 
tions (31) and (33). Such general solutions can be expressed by a 
standard procedure, and details of such cases will be reported in the 
future. 

The remaining part of this section will be devoted to a special 
case in which the absorptivity of the gas is virtually zero. This is 
the case of practical interest when the gas is optically thin. 

By putting w to be equal to unity, one finds that the solution 
satisfying equation (55) and the boundary conditions, equations (31) 
and (33), is given by: 
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( 57 ) 



4 n a r„ 


3 3 r" 



and 


16 ir a n r. 


Qr = 


3 6 



(58) 


The above expression, equation (58), can be substituted into equa- 
tions (44), (47), and (48) to predict the transfer number enhancement 
factors, e^, Cy^, and C^, as well as the radiation factor, equation (45). 
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4. 


PERFORMANCE ANALYSIS OF SELECTED LIQUID ROCKET ENGINES 


A numerical analysis for the SSME and VTE has been conducted to 
compare the performance characteristics of high and low enthalpy engines 
and to assess the need of regeneration of the heat transported to the 
chamber wall of each engine. 

In order to achieve the objectives, the sensitivity of each engine 
on the following factors is considered: (1) engine size, (2) Nusselt 

number, (3) emissivity, (4) regenerative efficiency, and (5) power level. 

The following two sets of assumptions are adopted as the reference 
for the performance comparison: 

(1) The overall vaporization efficiencies of bipropellants and the 
overall combustion efficiency, in the absence of radiation, are given by 


VTE : 

n VF = °- 94 

^VO 0.95 

and 

o 

ii 

o 

VO 

N5 

SSME : 

^VF - 

n vo 0-9 

and 

J=S 

n 

ii 

o 

.92 


(2) The overall vaporization efficiencies for fuel and oxidizer 
droplets and the global combustion efficiency in the liquid rocket 
engine combustion chamber are given by 


_ c 

n VFR ihyp n VF VFR n VF 


n VOR 


*VOR 

”Vo nvo 


C VOR n V0 


(59) 

jr 

\ 

(60) 


where Cy^ and C VQR are obtained from equations (47) and (48), respec- 
tively. These relationships are consistent with the classical criterion 
of liquid rocket engine design advocated by Priem (1958). For example, 

C„ ’ s are determined by 
VK 


c 


VR 


1 + (l-n R ) [C p (T m -T b )/L]) 


"‘1 

' i ■ 

■ ["it C P<W /L ] 

In 

. c p <v t b>1 


. L J 


(61) 
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The mean gas temperature, appearing in equation (61), is approximated by 



T + 5 T ) 


(62) 


in which and 5 are the weighing factors, respectively, depending on 

the structure of bipropellant combustion (Chiu, 1986a, 1986b). For 
example, large group flame emission is largely responsible for the heat 
transfer to the wall and to the upstream regions to enhance droplet 
vaporization. Envelope flames of individual droplets and droplet clusters 
serve to self-accelerate the gasification process. In the latter cases, 
the numerical value of of equation (62) is close to 2 and is zero. 

Vaporizing droplets or clouds near the wall will have 5 values close to 

2 and close to zero. The nature of the distributions of the group 

flame and drop or cluster bound flames strongly depends on the injection 
characteristics, droplet mixing pattern, and various spray-turbulent 
interaction processes. In view of this cqmplex mechanism and the lack 
of experimental data or numerical simulation at the present time, the 
average values of and will be taken to be unity, i.e.. 



T ) 
w 


4.1 Results and Discussion 

The scaling laws and gasification rate models are used to develop 
a numerical algorithm (see Appendix A), to aid in the systematic numerical 
prediction of the chamber temperature as a function of chamber design 
and operating conditions. Basic engine data are listed in Table 1. 

4.1.1 Effects of Chamber Geometry 

Since the radiative heat loss is approximately linearly proportional 
to the chamber volume, whereas the conductive heat loss increases linearly 
with respect to the chamber surface area, the characteristics of chamber 
temperature loss with respect to the geometry are investigated. Figure 
4 shows that the increase in the effective length and diameter of the VTE 
combustion chamber reduces the chamber temperature appreciably. Note 

that the maximum loss occurs at the highest gas emissivity, i.e., e = 1.0 

§ 

and q = 0. The maximum loss in a prototype engine is approximately 
100°R. However, when the effective length is increased by 10 times the 
real combustion value, i.e., = 3.33 ft, the maximum temperature 
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Fig. 4. Radiation Effect on the Performance of the VTE 
and Additional Hypothetical Sizes 
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reduction is approximately 600 °R. The ratio of the radiative loss to 
the conductive heat loss is approximately 2 to 1 . 

In contrast to the VTE, the chamber temperature reduction in SSME 
is negligibly small. For example, the maximum chamber temperature reduc- 
tion in a prototype engine, with the emissivity being equal to unity 
and a vanishing regenerative efficiency, amounts to a few degrees. This 
is in sharp contrast to VTE with a temperature reduction of 100°R as 
previously described. Figure 5 shows that the maximum temperature drop 
of 150°R could occur in a hypothetical SSME engine that has a 200 times 
larger chamber volume than the prototype. 

These significantly different chamber temperature loss characteris- 
tics of the VTE and SSME are attributed to the low and high enthalpy 
flow characteristics of each engine. 

Figure 6 shows the ranges of the variations of y and K for two types 
of engines, the VTE and the SSME. Note that the variations in the values 

-5 -3 

of y for the SSME range from 10 to 10 , and K remains nearly constant 
as the engine size increases. Note that the values of y in the VTE are 
approximately 100 times greater than that of the SSME and a pronounced 
variation in K value occurs when the engine size increases. The latter 
factor has a small positive effect for minimizing the chamber tempera- 
ture reduction. 

The results shown above and in the remaining part of this section 
suggest that the low-enthalpy engines could have a relatively large 
temperature drop due to radiative heat loss. The regenerative scheme 
will be effective in particular when such engines are designed or being 
operated in marginal conditions such as lower power level flight. 

4.1.2 Effects of Nusselt Number 

The relative importance of the conductive-convective heat loss to 
that of radiation is assessed by comparing the chamber temperature 
variations in a range of Nusselt number extending from an order of 
3 

unity to 10 . Figure 7 shows that the variation in T for the VTE is 

2 ^ 

relatively small in the range of 1 < N^ < 10 . A significant conductive- 

convective loss occurs when N^ exceeds a critical value of approximately 

300. The presence of such a critical Nusselt number suggests that the 
low-enthalpy engines, such as the VTE, may be further classified into the 
"radiation loss dominated engine” when N^ is smaller than the critical 

value and "conduction-convection loss dominated engine" if is greater 

than the critical value. 



Fig. 5. Radiation Effects on the Performance of the SSME 
and Additional Hypothetical Sizes 
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Fig. 7. Effect of Nusselt Number and Vaporization Efficiency on 

VTE Chamber Temperature 


In contrast to the VTE, the chamber temperature losses estimated for 
the SSME in the same range of Nusselt numbers are found to be negligibly 
small (Fig. 8). 

4.1.3 Effects of Radiation- Enhanced Combustion and Vaporization Rates 

The effects of radiation of the propellant gasification enhancement 
in the SSME is shown in Figure 9. Note that an almost 150°R gain in the 
chamber temperature is achieved in an adiabatic flame temperature case 

”4 _3 

when ri R Q is increased from 5 x 10 to 2 x 10 . The gain in chamber 

temperature at various degrees of regeneration is approximately 150°R 
for the prescribed increase in n Figures 7 and 8 also show the general 

feature of the chamber temperature gain in the range of the Nusselt numbers 
considered . 

This general temperature gain characteristic is attributed to the 
fact that the vaporization enhancement rate is approximately propor- 
tional to the enthalpy flows and, thus, constitutes a significant frac- 
tion of energy addition into the combustion process. 

Thus, for both high- and low-enthalpy engines, the stretegy is to 
use as much radiation energy as possible to enhance the gasification 
process and the associated combustion efficiency while minimizing the 
radiative transfer loss to the wall. 

4.1.4 Effects of Engine Power Level 

Since the throttling of the engine power output is accomplished, in 
general, by the reduction in total enthalpy flow, the lower power level 
operation will exhibit low enthalpy characteristics; i.e., high y engine 
behavior. Figure 10 shows the VTE chamber temperature variations in the 
ranges of power level and gas emissivity between 0 to 1.0. The effects 
of the Nusselt number, the power level and gas emissivity on the chamber 
temperature are shown in Figure 11. 

The results once again suggest that a dramatic reduction in chamber 
temperature, on the order of a few to several hundred degrees, is 
expected in low enthalpy engines. 
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SSME Chamber 
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Effect of Radiation on VTE Chamber Temperature at Selected 
Nusselt Numbers and at Various Power Levels 
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5. 


CONCLUSION 


The scaling law coupled with the radiation-enhanced propellant combus- 
tion model and associated computational algorithms have been developed 
to assess the impacts of radiation on the chamber temperature of liquid 
rocket engines. Two major impacts due to radiation are (1) the reduc- 
tion in the temperature due to unrecovered radiative conductive and 
convective heat transfer to the chamber wall with a radiative amount of 
approximately 70% of the total loss when the Nusselt number is smaller 
than a few hundred, and with increasing dominance by a conductive- 
convective loss at larger Nusselt numbers, and (2) the gain in the cham- 
ber temperature due to the radiation-enhanced gasification rate that 
tends to increase the combustion efficiency by a few percent over the 
non-radiating combustion environment. 

Based on the characteristics of radiative loss and gain, as pre- 
sented in the numerical illustration, liquid rocket engines can be 
classified into high-enthalpy engines, such as the SSME; and low-enthalpy 
engines, such as the VTE. The performance of the latter type of engines 
can be significantly impaired by radiative loss at all power levels with 
high emissivity due to heteropolar gaseous molecules such as H^O, CO^, 

CO, NH^, and NO^. A regenerative system with an efficiency higher than 

50% would be required to maintain an acceptable engine performance. 

Although the overall combustion efficiency model presented here is 
crude, and uncertainties of the numerical values adopted may affect the 
accuracy, the positive impact of radiation-enhanced droplet vaporization 
appears to be a fundamental consideration in the selection of the pro- 
pellant and the design of the injector. Much experimental data and 
extensive modeling, especially with CFD, are required to support the 
preliminary discoveries. 

The overall assessment of the general problem of radiation phenomena 
in liquid rocket engines clearly indicates the need for research in the 
following specific problem areas. 

(1) Droplet-radiation Interaction — The problems related with the 
radiative effects on the droplet ignition, combustion, extinction, as 
well as droplet saturation, condensation of cryogenic and non-cryogenic 
propellants, are of basic importance in radiation- coup led spray combus- 
tion. The effects of the collective interaction phenomena can signifi- 
cantly alter the radiation beam penetration, scattering, absorption, and 
emission. Droplets near walls, in envelope flames, and recirculation 
zones are subject to specific radiation view factors which can promote 
or impair the gasification rates. 
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(2) Turbulence-Radiation Interaction — High temperature, high 
intensity turbulence in the combustion zone provokes turbulence-radiation 
interaction when the radiative cooling takes part in viscous dissipation 
and heating. In particular, when the wavelength of the radiation is 
comparable to the scale of the turbulence, the intercoupling could 
contribute strongly to the radiative decay of turbulence. Many of these 
effects are anticipated to depend on the molecular properties, the 
chamber pressure, and temperature. Heterogeneous mixtures, soot par- 
ticles and droplets which have high emissivity and absorptivity are 

also important factors in the local turbulence-radiation interaction. 

(3) Flame-Radiation Interaction — Thermal emission and chemi- 
luminescense in a flame zone are the cause of the radiative cooling in 
optically thin flames. The cooling impairs the flame stability charac- 
teristics. Many of these phenomena depend on the molecular and thermo- 
chemical properties of the gas mixture in the preheating and combustion 
zones. In assessing the interaction between the flame and radiation, 
one needs to differentiate premixed versus diffusion/mixing controlled 
flames because the heat release rate and the cooling rate in each type 
of flame are characteristically different. This suggests that the 

flame extinction characteristics in the two types of flames are different. 
Additionally, the presence of soot particles in the flame environment 
affects the emissivity and absorptivity of the mixture. Systematic 
studies of these fundamental processes under typical rocket engine 
operating conditions are still lacking at present. 

The prediction of combustion processes in a radiative environment 
requires an advanced spray combustion CFD code that is coupled with a 
radiative transfer subcode. The latter accounts for the non-equilibrium 
radiation transport processes. In addition, combustion submodels 
including radiative droplet laws, atomization processes, turbulence 
models, gas phase combustion, as well as emissivity, absorptivity models 
in super or subcritical states under radiation effects must be supple- 
mented. Finally, new computational techniques based on finite difference 
schemes for conventional CFD and special computational techniques such 
as multi-flux schemes must be effectively incorporated. In addition, 
the initial and boundary conditions must be supplied to formulate a 
well-posed computational problem. 
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APPENDIX A 


PROGRAM LISTING 


ORIGINAL PAGE IS 
OF POOR QUALITY 


PAUL CHIU'S RADIATION COMPUTER MODEL 

DIMENSION ETAREG( 11 ) ,GRAD( 10 ) , INCL( 2 ) , ETARF (10), ETARO ( 10 ) 

DATA FTLBTU/778. 16/, GC/32. 17 4/, KG, LE,MF,NO/4*l/, PI/3. 141592653/, 
1 IP,IRAD/0,0/,MAX,MAXR/2*50/,TOL/0.0020/,INCL/'EXCL' , ' INCL ' / 
REAL ITLIM, KCOND , KP , L ENGTH , LF , LO, MDOTF , MDOTO , NUS S EL 

NAMELIST /PAUL/ CPFI , CPOI , CPPROD , DIAM , ETAC , ETAREG , ETARF , ETARO , 


ANK 8/88 
ANK 8/88 
ANK 8/88 
ANK 7/88 


/PAUL/ 

ETAVF , ETAVO , ETAWAL , GRAD , IP , I RAD , ITLIM , KCOND , KG , LE , /PAUL/ 
LENGTH , LF , LO , MAX , MAXR , MDOTF , MDOTO , MF , NO , NUSSEL , /PAUL/ 
QCOMB , SIGMA , TBF , TBO , TFI , TOI , TOL , TWALL , ALPHA , PHI , QE /PAUL/ 


READ (5, PAUL) ANK 7/88 

PHI - MDOTO/MDOT F ANK 7/88 

PHI 1 = 1.0 + PHI ANK 7/88 

QE - ETAC*PHI 1 *QCOMB - ETAVF *LF - ETAVO*LO*PHI ANK 7/88 

ALPHA = KCOND*LENGTH*NUSSEL*PI/(2.0*CPPROD*FTLBTU*GC*MDOTF*PHIl) ANK 7/88 

WRITE (6,2) INCL( IRAD+1 ) ANK 8/88 

2 FORMAT (//2 OX, 'PROPELLANT VAPORIZATION AFFECTED BY REGENERATED HEAANK 8/88 


IT (ETAREG) ',A4,'UDING RADIATION'//) 

ANK 

8/88 

WRITE (9,2) INCL( IRAD+1) 

ANK 

8/88 

WRITE (6, PAUL) 

ANK 

7/88 

WRITE (9, PAUL) 

ANK 

8/88 

DO 40 M - 1,MF 

ANK 

8/88 

DO 40 N - 1 , NO 

ANK 

8/88 

DO 40 K ** 1 , KG 

ANK 

8/88 

WRITE (9,6) 

ANK 

8/88 

FORMAT (/) 

ANK 

8/88 

DO 40 L - 1 , LE 

ANK 

8/88 

J » 0 

ANK 

8/88 

TCI - 0.0 

ANK 

8/88 


BETA - DIAM*GRAD(K)*LENGTH*PI*QE**3*SIGMA/( 3600 . 0*CPPROD**4*MDOTF ANK 8/88 
L *PHIl**4 ) ANK 7/88 
ALPHAP - ALPHA* (1.0 - ETAREG ( L ) ) ANK 8/88 
BETAP - BETA* (1.0 - ETAREG ( L ) ) ANK 8/88 


BETAP/( ALPHAP + 1.0) 

(CPFI + CPOI*PHI )/PHIl 

( CPFI *TFI + CPOI * PHI *TOI)/(CPFI + CPOI*PHI) 

CPROPI *PHIl*TPROPI/QE 
CPPROD* PHI 1 *TWALL/QE 

+ ALPHAP*THETAW + BETAP*THETAW**4 + THETAI )/( ALPHAP+1 . ) 


10 


I .GT. MAX) GO TO 30 


20 


30 


GAMMAP 
CPROPI 
TPROPI 
THETAI 
THETAW 
KP - (1 
QC - KP 
I - 1 

REST - GAMMAP*QC**4 + QC - KP 
IF ( ABS ( REST ) .LT. ITLIM .OR. 

IF ( I .GE. 2) GO TO 20 
QC1 - QC 
RESTl « REST 
1*2 

QC - 1.2 0*QC 
GO TO 10 
QC2 - QC1 
REST2 - RESTl 
QCl - QC 
RESTl - REST 

IF ( ABS ( REST2 - RESTl) . LE . 0 . 0 ) GO TO 30 
QC - QCl - RESTl* (QC2 - QCl)/(REST2 - RESTl) 
I - I + 1 
GO TO 10 

TC - QC*QE/( CPPROD*PHIl ) 

IF (IRAD .EQ. 0) GO TO 35 
RADIATION IS TO BE INCLUDED 
TC2 - TCI 
TCI - TC 


ANK 8/88 
ANK 7/88 
ANK 7/88 
ANK 7/88 
ANK 7/88 
ANK 8/88 
ANK 7/88 
ANK 7/88 
ANK 7/88 
ANK 7/88 
ANK 7/88 
ANK 7/88 
ANK 7/88 
ANK 7/88 
ANK 7/88 
ANK 7/88 
ANK 7/88 
ANK 7/88 
ANK 7/88 
ANK 7/88 
ANK 7/88 
ANK 7/88 
ANK 7/88 
ANK 7/88 
ANK 7/88 
ANK 8/88 


ANK 

ANK 


8/88 

8/88 
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I 


35 


40 

4 


PROGRAM LISTING (Continued) 



J - J + 1 

ANK 

8/88 

TRATIO - ABS( TC2 - TC1)/TC1 

ANK 

8/88 

IF (TRATIO .LE. TOL .OR. J .GT. MAXR ) GO TO 35 

ANK 

8/88 

TCR - (TC + ETAWAL* TWALL )/2.0 

ANK 

8/88 

AF - CPPROD* (TCR - TBF )/LF 

ANK 

8/88 

AO - CPPROD* (TCR - TBO)/LO 

ANK 

8/88 


CRADF - LOG( ( 1 . 0+ ( 1 . 0-ETARF ( M ) ) *AF )/( 1 . 0-ETARF ( M) *AF ) )/LOG( 1 . 0+AF )ANK 
CRADF -1.0 FOR 

CRADO - LOG ( ( 1 . 0+ ( 1 . O-ETARO ( N ) ) *AO)/( 1 . 0-ETARO( N ) *AO) )/LOG( 1 . 0+AO)ANK 
ETACR - (CFADF*CRADO)*ETAC ANK 

ETAVFR - CFADF*ETAVF ANK 

ETAVOR - CFADO*ETAVO ANK 

QE - ETACR* PHI l*QCOMB - ETAVFR*LF - ETAVOR*LO*PHI ANK 

GO TO 5 ANK 

IF (IP .GT. 0) GO TO 40 ANK 

WRITE (6,1) I , K, L,CPROPI , ETAREG( L ) ,GRAD( K) , KP,QC, TC,TPROPI , ANK 

1 ALPHAP , BETAP , GAMMAP , REST , THETAI , THETAW ANK 

FORMAT (/' I -',14,' K -',13,' L -',13,' CPROPI -',F9.6, ANK 

1 ' ETAREG -' ,F5.3, ' GRAD -' , F5 . 3 , ' KP-',F9.6,' QC-',F9.6, ANK 

2 ' TC - ' , F7 . 1 , ' TPROPI -',F7.1/' ALPHAP- ', 1PE13 . 6 , ' BETAP-', ANK 

3 E13.6,' GAMMAP-' ,E13. 6, ' REST- ' , E13 . 6 , ' THETAI- E13 . 6 , ANK 

4 ' THETAW-' , El 3. 6/) ANK 

WRITE (6,3) J , M , N , CRADF , CRADO , ETACR , ETAVFR , ETAVOR , QE , TCR , ETARF (M) , ANK 

1 ETARO(N) ,AF,AO,TRATIO ANK 

FORMAT (' J-' , 13, ' M- ' , 1 3 , ' N-',I3,' CRADF- ' , F8 . 5 , ' CRADO-', ANK 

1 F8 . 5, ' ETACR-' ,F8.5, ' ETAVFR-' , F8 . 5 , ' ETAVOR- ', F8 . 5 , ANK 

2 ' QE- ' , Fl 0 . 3 , ' TCR- ' , F8 . 2/' ETARF(M) -',F9.6,' ETARO(N) -' , ANK 

3F9.6,' AF - ' , F10 . 3 , ' AO-',F10.3,' TRATIO - ' , El3 . 6/) ANK 

WRITE( 9 , 4 ) L , ETAREG ( L ) , TC , GRAD ( K ) , ETARF ( M ) , ETARO ( N ) , TCR , QC , QE ANK 


FORMAT ( /' L -' ,13, ' 

1 ' ETARF - ' , F8 . 6 , ' 

2 ' QE -' , F10 . 3 ) 

END 


ETAREG 
ETARO - 


,F6.4, 
,F8.6, ' 


TC 

TCR 


,F8.2, 
» F8 . 2 , 


GRAD 
QC - 


, F6 . 4 , ANK 
,F9.6, ANK 
ANK 
ANK 


8/88 

SSME 

8/88 

8/88 

8/88 

8/88 

8/88 

8/88 

8/88 

8/88 

8/88 

7/88 

8/88 

8/88 

7/88 

7/88 

8/88 

8/88 

8/88 

8/88 

8/88 

8/88 

8/88 

8/88 

8/88 

8/88 

7/88 


INPUT DATA 


$PAUL 

CPFI - 0.70, CPOI - 0.3780, CPPROD - 1.650, DIAM - 0.20412, ETAC - 0.920, 
ETAREG - 1.,.9,.8,.7,.6,.5,.4,.3,.2,.1,.0, ETAVF - .94, ETAVO - .95, IRAD - 1 
GRAD - 1.0, 0.8, 0.6, 0.4, 0.2, 0.0, ITLIM - 1.0E-10, MF - 1, NO - 1, 

KCOND - 2.0, KG - 6, LE - 11, LENGTH - 0.33333, LF - 377.0, LO - 178.20, 

MDOTF - 0.1415728, MDOTO - 0.2345336, NUSSEL - 300.0, QCOMB - 8910.0, 

SIGMA - 1.7120E-9, TFI - 460.0, TOI - 460.0, TWALL - 2500.0, TOL - 0.0020, 
ETAWAL - 1.0, ETARF - .008,. 01,, ETARO - .0005,. 001, TBF - 649.0, 

TBO - 530.0, 

$END 
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